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We analyze the full counting statistics of photons emitted by a double quantum dot (DQD) to a 
high-quality microwave transmission line due to the dipole coupling. We show that at the resonant 
condition between the energy splitting of the DQD and the photon energy in the transmission line, 
photon statistics exhibits both a sub-Poissonian distribution and antibunching. In the ideal case, 
when the system decoherence stems only from photodetection, the photon noise is reduced below 
one-half of the noise for the Poisson distribution. The photon distribution remains sub-Poissonian 
even at moderate decoherence in the DQD. 

PACS numbers: 73.23.-b, 73.63.Kv, 42.50.Ar 



Introduction. The statistics of photons emitted by 
an electric current depends on an electron state of a con- 
ductor. If electric current were classical, photon field 
would be in a coherent state^ with the Poisson noise. A 
quantum electron system with strong inelastic processes 
is characterized by the thermal distribution and produce 
black-body radiation with super-Poissonian statistics of 
emitted photons. However, if the electron distribution is 
far from equilibrium, the photon counting statistics may 
become sub-Poissonian.-'' ' 

Several experimental approaches have recently been 
developed to study statistics of photons in the GHz 
frequency range. Experiments''' measure the photon 
statistics in a steady state of high quality resonator and 
distinguish between the thermal source and a coherent 
drive. Photon noise of a quantum point contact at finite 
bias was also investigated using the amplifier technique.*' 
An alternative approach to study photon statistics uti- 
lizes photon counter.''^ An individual photon counter can 
provide information about statistics of emitted photons, 
while a system with two counters can be used to measure 
the second order intensity correlation function g^^ (r) , 
see e.g. Ref. 7, which represents the correlations in ob- 
serving two photons at two moments separated by time 
r. 

In this paper we analyze the statistics of photon radia- 
tion from a double quantum dot coupled to a high-quality 
microwave transmission line, a system that was recently 
experimentally studied by several groups.''"'^ In the ab- 
sence of decoherence processes in the DQD, the transition 
from its excited state with energy Eg to the ground state 
with energy Eg occur only via emission of a photon to 
the transmission line. In this case, dependences of the 
current and the average photon number on energy dif- 
ference Ee — Eg exhibit resonant peaks at photon energy 
hujo = Ee — Eg. At the same condition, we find the sup- 
pression of noise in charge transfer and photon emission, 
indicating on sub-Poisson statistics. 

We also evaluate the second order correlation function 
■9ph (''')■ ^^'^ ^^^^ ^^^^^ function is reduced for time 
intervals shorter than the typical time of an electron pas- 
sage through a DQD, thus reflecting an electronic origin 
of photon emission in the resonator. We also calculate 
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FIG. 1. (Color online) (a) An illustration of a DQD and 
a transmission line coupled to a photon counter (PC), (b) 
In the DQD, electrons are confined to the left (L) and right 
(R) dots by barrier gates BL, BM, and BR that also control 
electron tunneling rates between the source, S, and the left 
dot, the left and right dots, and the right dot and the drain, D, 
respectively. Electrostatic energies of two quantum dots are 
defined by the plunger gates, PL and PR, and the PL gate 
is also connected to an antinode of the transmutation line, 
(c) Electronic states of the DQD are presented in both the 
eigenstate basis (solid lines) and the left-right basis (dashed 
lines). Tunneling from the left lead to the ground/excited 
state with rate Fi^g/e, and from the ground/excited state to 
the right lead with rate Fg/^ are illustrated by arrows. 



the distribution function for photon counting, which dif- 
fers significantly from the Poisson distribution. Finally, 
we analyze the effect of energy and phase relaxation of 
electron states in the DQD. We find that while pure de- 
phasing has a weak influence on the photon statistics, in- 
elastic energy relaxation changes photon counting statis- 
tics towards super-Poissonian. 

Model. We study statistical properties of charge trans- 
fer through the DQD and photon emission to the trans- 
mission line. The Hamiltonian for a system of coupled 
DQD and the transmission line, shown in Fig. 1, is pre- 
sented as a combination of three terms, H = iJoQD + 
^ph + ffint- The non-interacting DQD near a triple point 
in its electrostatic stability diagram ' - is represented by 
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in the basis of electron states in the left, \L), and 
right, \R), quantum dots with electrostatic level bias 
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FIG. 2. (Color online) Dependences of (a) the electric current 
/ through the DQD, (b) the average photon number A'^ in the 
transmission line, (c) the electron Fano factor F^i and (d) the 
Fano factor Fph for emitted photons are shown as functions 
of the electrostatic bias he near the resonance at teo with 
£0 ~ \/ uj'i — iT^ ■ Solid (dashed) lines represent the case for 
7r = (7r ~ 2r). Other system parameters are Fi^^ ~ F, 
Luo = 800F, T = 200r, go = 5r, k = 2F and 7^ = 0. Dotted 
line in panel (a) refers to elastic electric current through a 
non-interacting DQD. 



e and the tunneling amplitude T; in this basis, = 
\R) (R\ - \L) {L\ and = \R) {L\ + \L) {R\. The term 
TJph = hujoa^a represents a noninteracting photon mode 
in the transmission line. The interaction between charge 
and photon degrees of freedom is described by the Jaynes- 
Cummings Hamiltonian Hint ~ %o (aV~ + acr+) 

In further calculations, we use the eigenstates of the 
DQD Hamiltonian, Eq.(l), namely the ground, \g), and 
excited, |e), states: 



|e) = cos(6l/2) \L) + sin(6i/2) \R) , 
\g) = - sin(6l/2) \L) + cos(6i/2) \R) 



(2) 



Here 9 = arctan(2T/e) characterizes the hybridization 
between states \L) and \R). The energy splitting be- 
tween the eigenstates hft — hy/e^ + 4T^ can be tuned 
via gate voltages. In the eigenstate basis, Eq. (2), the 
full Hamiltonian is reduced to 



H 



+ kbjQo} a + hg[a} a + acr"'"), (3) 



where g = go sin 9 is the effective electron-photon cou- 
pling constant, and the matrices a~ = \g) (e|, cr"*" = 
|e) {g\ and = \e) (e| — \g) {g\ are defined in the eigen- 
state basis. 

We analyze the behavior of the system with Hamilto- 
nian Eq. (3) in the presence of decoherence in electron 
and photon degrees of freedom and tunneling of electrons 
between the DQD and the leads by employing the Born- 
Markov master equation for the full density matrix 



Here, the first term on the r.h.s. of Eq. (4) describes 
the unitary evolution of the system and the second term 
accounts for the decoherence and tunneling processes in 
the DQD and the transmission line: 



VtotP=K'D{a)p + Jr'D{<T )p + ^X>(crz)/0 



(5) 



where T){x)p — {2x px^ x p ^ px^ /2 is the Lind- 
blad superoperator. Detection of a photon in the trans- 
mission line by a photodetector with rate k is represented 
by the first term, KD{a)p and the electron relaxation from 
the excited state |e) to the ground state \g) with rate 7^ 
is represented by ^rT^{(j~)p. Similarly the third term 
describes the dephasing of the DQD with rate 7^. 

The last two terms in Eq. (5) account for the processes 
of loading state \L) from the source with tunneling rate 
F; and unloading state \R) to the drain with tunneling 
rate F^, see FIG. 1. The tunneling superoperators are 
expressed in terms of Cr — |0) and c\ = \L) 

For this model, one can evaluate both charge trans- 
fer and photon emission statistics. The charge transfer 
statistics has been discussed in the literature for similar 
systems^''' ^'"^^ and we present main relations for charge 
transfer statistics through a DQD in the Supplementary 
Materials. Below we provide equations for photon statis- 
tics. We assume that photon detection is performed by 
an ideal detector so that all emitted photons arc even- 
tually captured by the detector. The latter assumption 
is justified for high quality superconducting transmission 
line which can be made practically lossless. 

The probability distribution Pn(^o) to count n photons 
during measurement time tg is related to the generating 
function (GF) C/(to,s) through a discrete Fourier trans- 
form with respect to counting field x G [0, 27r]: 



e(to,e-^'^)e^'^" J. 
27r 



(6) 



The GF G{t,s) = Tr{p(t, s)} is expressed in terms of 
the generalized density matrix p{t, s) = e^^*)*pst, which 
is the solution of the master equation with parametric 
dependence on the counting field s = e~'-^: 

= M{s)p{t, s), M{s) =C+(s- l)J{a). (7) 

Here J{a)p = napa^ is the jump operator that describes 
the photon detection with transition rate k and pst is the 
steady state of the density matrix defined by Eq. (4) with 
I^Pst = 0. 

To quantify the width of the distribution P„(io)j we 
evaluate the photon noise Fano factor 



ph 



(8) 



Lp 



[H, p] + Ptotp. 



(4) 



The Fano factor equals 1 for Poissonian processes, while 
for sub-(super-) Poissonian processes, the Fano factor is 
bclow(abovc) Ir 
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FIG. 3. (Color online) (a) The second order correlation func- 
tion (;pij(2)(r) for photons as a function of time r at resonant 
condition — ujo for different values of energy, 7^, and phase, 
70, relaxation rates, (b) The photon Fano factor Fph as a 
function of relaxation rate 7 in three cases 7,. = 7,^ = 7 (solid 
line), 7r = 7 and 7,/, = (dotted line), 7,/) = 7 and 7r = 
(dashed line), (c) Fph as a function of photon detection rate 
K shows a flat behavior for k > F. Other system parameters 
in both panels are F;,,. = F, wo = 800F, T = 200F, go = 5F, 
K = 2F and 7^ = O[for (a) and (b)]. 

For a lossless transmission line, the average photon 
count (n) = uNto during time to is determined by 
the photon number N = (a^a)st — Tr{psta^a} in the 
transmission line and the photon detection rate k.-~ If 
the measurement time exceeds the characteristic system 
memory time, the Fano factor is <o~ independent and can 
be evaluated as a time integral:" ' 

Fph = l + 2«;7V^ (4h(r)-l)dr. (9) 

Here, the second order correlation hmction ffpjj (t) is 
given by' 




(10) 

This function characterizes the joint probability to ob- 
serve two photons separated by time interval r and its 
asymptotic value is 1 at long times. In particular, photon 
antibunching (bunching) occurs when 5pjj (0) is smaller 
(greater) than 1. 

We discuss alternative methods of calculating statisti- 
cal properties of photon emission in Sec. S.I. of the Sup- 
plementary Materials. We note that one can also apply 
similar formalism to evaluate cross-correlation functions 
of electron charge transfer and photon emission, which 
will be discussed elsewhere. Experimental observation of 
such cross-correlations is a challenging task, but can be 
achieved by combining charge sensing measurements- ' 
and photon detection.'" 



Results. Below we consider a DQD with equal tunnel- 
ing rates through left and right contacts, Ti^r = T, and 
the interdot tunneling amplitude T = 200r. We take 
Wo = 800F, the photon relaxation rate k ^ 2T and the 
electron-photon bare coupling go = 5F. 

In FIG. 2 we present dependence of (a) electric cur- 
rent /, (b) the average number of photons in the trans- 
mission line, N = {a^a)st, (c) the Fano factor of charge 
current, F^i, and (d) the Fano factor for photon flux, 
Fph, Eq.(8). Solid lines in FIG. 2 are evaluated for an 
ideal quantum dot when the inelastic relaxation from 
excited to the ground state is absent. In this case, 
the amplitudes of electric current and the photon flux 
have a well pronounced peak at the resonant condition 
f2 = Wo, while away from the resonance, photon produc- 
tion is suppressed, iV — >■ 0, and the current approaches 
lo = eT^r/(e^+3T^) for elastic electron transfer through 
a non-interacting DQD,-'' shown by a dotted line in 
FIG. 2a. Fano factors for both electric current and pho- 
ton flux are reduced below 1/2, indicating sub-Poissonian 
statsitics with strong suppression of charge and photon 
noise at the resonance. 

Moderate inelastic relaxation changes the above pic- 
ture. Because inelastic relaxation facilitates electron 
transfer through the DQD, the electric current exceeds Iq 
even away from the resonance, \n — wpj ^ F, and only 
a weak enhancement of the current occurs a,t Q = ujq. 
The electron Fano factor is reduced below unity in the 
presence of energy relaxation and the resonant electron 
transfer with photon emission does not significantly af- 
fects Fci- We still observe a resonant emission of photons, 
see dashed line in FIG. 2b, but the statistics of emitted 
photons is closer to Poissonian. 

We study the properties of the second order correlation 
(2) 

function (''")' FIG. 3a. The thick solid line shows 

(2) 

9ph (''') ^'^^ ideal DQD without inelastic relaxation and 
dephasing, 71. = 7^ = 0. Probability to observe two 

(2) 

photons simultaneously is reduced, .gp^ (0) < 1, indicat- 
ing photon antibunching. As r becomes longer than the 

typical electron transfer time through the DQD, ^ l/F, 

(2) 

function g^^^ (r) increases and eventually approaches its 
(2) 

asymptote, ffp^ (t — > 00) = 1. The integral in Eq. (9) 

(2) 

with such gfpj^ (r) is negative and the photon noise Fano 
factor Fph < 1. 

In the presence of inelastic relaxation in the DQD, 
Spti (*-*) increases and (r) reaches its long-time asymp- 
totic value 1 at shorter time scale. This behavior of 
(2) 

•9ph (''') results in larger value of Fp^. Pure dephasing, 
70, only smoothens <?ph (''')j and has little effect on the 
integral in Eq. (9) and, consequently, on i^ph- 

To investigate dependence of -Fph on inelastic, 7,-, and 
dephasing, 7^, rates, we plot -Fph(7r,70) as a function of 
7r. and 7^ in FIG. 3b. Pure dephasing, 7^ = 7 and 7^ = 
0, has extremely weak effect on J^ph even for large values 
of 70, while the inelastic relaxation, 7^ = 7 and 7^ ~ 0, 
quickly recovers i^ph to its Poissonian value, i^ph = 1. 
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FIG. 4. (Color online) Probabilities P„ to have n emitted 
photons during time to at resonance Q = ujo for (a) an ideal 
DQD without decoherence of electronic states, 7r = 7,^ = 
and to = 75/r; (b) a DQD with inelastic relaxation 7^ — 2r, 
70 = and to = 145/r. Other system parameters in both 
panels are T;,^ = T, a;o = SOOF, T = 200r, go = 5r and k = 
2r. A thin curve in both panels represents the corresponding 
Poisson distribution Pn^' = e~^n^ ln\ with n = nP^ 
equal to the average number of emitted photons. 



Again, additional pure dephasing, 7^ = 7r = 7, makes 
no significant corrections to i^ph- The weak dependence 
of i^ph(7r,70) on 70 is consistent with weak dependence 
of gpjj (r) on 7^, as discussed above. 

We note that another Fano factor can be introduced to 
characterize the variance of the photon number N ^ a^a 
in the steady state, F^^f = ((TV^^.t - (iV)2j/(iV),t. ^^The 
quantity F^^ = 1 + (iV)st(.gph (0) - 1) is given by instan- 

taneous value .^pj^ (r = 0) in a steady state and does not 
represent time evolution of the photon field. Therefore, 

fst) 

Fpjj < 1 is an indicator of the photon antibunching. On 
the other hand, photon antibunching or bunching does 
not imply that the full counting statistics is sub- or supcr- 
Poissonian, as the type of statistics of emitted photons 
is determined by the time integral of the second order 
correlation function for photons, Eq. (9)."' 

Next, we study dependence of photon noise Fano fac- 
tor on the photon detection rate. As the photon de- 
cay rate decreases, the average photon number in the 
transmission line increases. At large A'^, photons already 
present in the resonator cause stimulated emission by the 
DQD '■ . For an ideal DQD without energy relaxation, 
7r = 0, the photon Fano factor grows fast for k < F, 
see FIG. 3c. On the other hand, if the energy relaxation 
in the DQD is significant, the photon Fano factor can 
exceed unity and the photon field exhibits properties of 
a thermal state. In this respect, inelastic processes in 
the DQD enhance photon noise. For strong coupling of 
the transmission line to the photodetector, k ^ 7rj 5, and 
consequently low photon number in the transmission line. 



the back action of the photon field on electrons is negli- 
gible. The adiabatic elimination method' can be applied 
to study photon statistics emitted by a quantum dot, see 
6.17. Ref. 29 and the Supplementary Materials. 

Finally, we present the distribution function of the 
photon counts n over time <o- For ideal DQD with 
7r = 70 = 0, we take — 75/r and obtain the dis- 
tribution Pn{to) from Eq. (6), shown in FIG. 4a. The 
average value of photon counts, (n) ~ J2n "^^n — 27.6 is 
consistent with {n) = NkIq with N = 0.184. The vari- 



ance of photon coimts, 



11.6, in turn gives 



the Fano factor F-ph = 0.42, which coincides with previ- 
ous calculation, see Fig. 2d. For comparison, we present 
the Poisson distribution with the same expectation value 
(n) by narrow black dots in FIG. 4a. Once we consider 
inelastic processes in the DQD, the distribution is closer 
to the Poisson distribution, as shown in FIG. 4b, where 
the vertical bars represent P„ for a system with 71. = 2r, 
to = 145/r, (n) ~ 27.6 {N = 0.095) and Fph = 0.76. 

Conclusions. We investigated statistics of photons 
emitted by a biased DQD coupled to a lossless trans- 
mission line. We calculated the time correlation function 
(2) 

dphi''') ^^^^ found that the probability to observe two 
photons simultaneously is smaller than that for longer 
times, indicating photon antibunching. We also calcu- 
lated photon counting statistics Pn{ta) of observing n 
photons during a fixed time interval to. We find that 
distribution P„(to) shows a sub-Poissonian statistics. 

In recent experiments, decoherence rates are compara- 
ble to the strength of the electron-photon coupling. For 
this reason, we investigated the effect on charge and pho- 
ton statistics of pure dephasing in the DQD and energy 
relaxation. We found that pure dephasing does not sig- 
nificantly modify the charge transfer or photon emission 
statistics, but the inelastic relaxation processes result in 
several drastic changes, see FIG. 2: (i) The electric cur- 
rent and its noise acquire strong background as the in- 
elastic processes facilitate the charge transfer throughout 
the DQD, and the peak in / and the deep in F^i at the 
resonant condition E^ — Eg = Hloq are flattened, (ii) The 
photon number N at the resonance is suppressed as the 
effective photon source is reduced due to additional chan- 
nels for e — >■ (7 transition via inelastic events. (Hi) Photon 
Fano factor i^ph as a function of the level spacing is flatten 
as well. In the presence of inelastic electron relaxation, 
the memory time of the DQD is reduced, that increases 
the photon correlation function ^pj^ at short time scales, 
see FIG. 3a, and brings photon distribution function to 
the Poissonian distribution, FIG. 4b. 
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Supplementary Materials for "Full Counting 
Statistics of Photons Emitted by Double 
Quantum Dot" 

S. I. PHOTON EMISSION STATISTICS 

S. I. 1. Quantum Jump Approach 

In this section wc utilize the quantum jump approach''" 
to an open quantum system. Tlie Liouvillian for the 
equation of the motion of the density matrix in Lindblad 
form can be decomposed as 

p(t) = /:p(i) = (/:o + ^)p(<), (sii) 

where we have defined the jump superopertor, J p = 
Kapa\ to describe the quantum jump that corresponds 
to detection of a photon, and Co governing the rest of 
the dynamics of the system. The full density matrix p{t) 
can be resolved in terms of its components p'"^ (t) repre- 
senting a process with n photons being counted by the 
photon detector during a time interval t: 

n 

The equation of motions for p^'^^ (t) is 

p^"\t)^CoP^"\t)+Jp^"~'^t). (S13) 

These equations of motion are coupled and therefore are 
hard to solve. It is more convenient to define the gener- 
alized density matrix 

p(i,s)=E,sV"Hi), (S14) 

n 

by introducing the counting variable for photons, s. The 
equations of motion for the generalized density matrix 
p{t, s) is obtained by multiplying Eq.(S13) by s" and tak- 
ing sum over n, 

Kt,s)=M{s)p{t,s), (S15) 

with 

Mis)=Co + sJ. (S16) 

For s = 1, Eq. (S15) reduces to the master equation 
p{t) = Mil)pit). The formal solution of Eq. (S15) is 

p(t,s)=e^(^)*p(0,5), (S17) 

where /5(0, s) = pst is the initial state. 

We introduce moment generating function 

g{t, s) = Tr {i5{t, s)} = Tr {e^(^)*p(0, s)} . (S18) 

This function permits one to calculate the full counting 
statistics, i.e. the probability distribution P„(t) of emis- 
sion of n photons during time interval t. Indeed, the n 



resolved density matrix allows us to obtain the full count- 
ing statistics of the system by taking trace of p^'^\t): 

P„(t) = Tr{p(")(i)}. (S19) 
Then, according to Eqs. (S14) and (S18), we identify 

= Es"P„(t). (S20) 

n 

The probability distribution P„ (t) is given by the inverse 
Fourier transform in parameter s ~ exp(ix): 

Pn{t) = J^^ e-'"^g{t, e'^)^. (S21) 

The factorial moments of n can be obtained by deriva- 
tives of g{t, s): 

-E^«WnV-). (S22) 

s=l n 1=0 

We note that Eq. (SI 7) is understood as a Dyson series, 
therefore the density matrix can be unravelled into a sum 
of contributions corresponding to n = 0, 1, 2, . . . photon 
detections according to 

p{s,t)=S{t,0)pst+y2 f dtn--- I dti (S23) 

where S{tit2) = exp [£o(^i - ^2)] ■ 

S. I. 2. Numerical Calculation of the Fano Factor 

In principle, the method descibed above can be used 
not only for analytical calculation of the Fano factor, but 
for numerical computation as well. However, as shown in 
Eq. (S22), evaluation of factorial moments involve deriva- 
tives of generating function Q{t,s) over s, which is not 
accurate and neat enough in practice for numerical cal- 
culations. In this subsection, we describe a numerical 
method more suitable for numerical evaluation of Fano 
factors that we used in the main text. 

Our aim is to calculate the photon counts of the system 
during measurement time to in the steady state, 

{5n{tof) = ]- r dh r dt2{{5J{h),5J{h)}) , 
^ Jo Jo 

(S24) 

where 8J{t) = J{t) - J{t) and J{t) nNit) = TrJ^p^t 
is the average photon count rate with n being the mea- 
surement rate of photons; {A, B\ stands for an anticom- 
mutator. We can take derivatives of Eq. (S23) over s, 
and find the correlation function of photon counts dur- 
ing measurement time to in the form (n(to)) — Jq° dtJ{t) 



S7 



and 



to 



dti / dt2 



S^Tvp{to,s) 



Jo Ss{ti)5s{t2) 
Ss{t) 

2 [" dh I ' dt2{J{ti)S{ti,t2)J{t2)) 
JQ Jo 



(S25) 



(n{to)). 



where TZ = QC is the inverse of the Liouvillian pro- 
jected out of the steady state. 

Equation (S28) is the key step in evaluating the photon 
noise spectrum. We also have to find the inverse of the 
Loiouvillian, TZ, and project the result out of the steady 
state. In practice, inverse of the Liouvillian matrix with 
large dimension is numerically unstable, but we can eval- 
uate the combination \ W)) = TZJ'\st)) determined by the 
following equation 



Note that we have implied ti > t2 in the second line and 
thus the term with ti < t2 should be added to symmetrize 
the expression with switching on time labels. Then, one 
can integrate Eq. (S25) with respect to ti and t2 



invito)) - / / dhdt2J^ 
Jo Jo 

(-to 



(S26) 



(n{to)) + 2 / dT{to - t) [Tr (J5(t) Jp^t) - J' 
Jo 

(nito)) + 2J2 f " dT{t„ - t) (.9(2) (r) - l) , 



and arrives at the celebrated Mandel's photon count- 
ing formula.- ' Taking into account that to is large com- 
pared to the characteristic memory time of the system, 
Eq. (S26) reduces to the expression for the photon Fano 
factor, independent from tg 



F = 



1 + 2J^ dTTr|,g(2)(r) - l| . (S27) 



Following Ref. 20, we introduce the "Dirac notation" 
in Liouvillian space for steady state |st)) = p^t and a dual 
vector ((e| = 1. The inner product defined in Liouvillian 
space is the trace over the "ket" in state "bra". For exam- 
ple, the inner product of the two former objects is given 
by ((e|st)) = Trpst = 1- It is then useful to define the 
projector V =1^^ ^ |st)) ((e| onto the steady state as well 
as its complement Q^l — V. Note that a useful property 
of V is CV = C\st)){{e\ = and VC = 0, and therefore 
C = Q£Q. The propagator S{t) in the second hue of 
Eq. (S26) can be decomposed as S{t) ^ V + QS{t)Q, 
and thus TriJVJpst) = ((e| J^|st)) ((e| J^|st)) = J\ We 
obtain 



F = 1 



dTjQSiT)QJ 



(S28) 



= 1- -TT{jQC-^QJpst} 
= 1- l((e\jnj \st)), 



C\W)) = CTZJ\st)) = QJ\st)) 

= J|st))-|st))((e|J|st)), 



(S29) 



where the second equality is obtained by the relation 
Cn = CQC-'^Q = £(1 - V)C-^Q = Q.' To this end, 
the solution \W)) in Eq. (S29) is equivalent to the in- 
verse TZ. In the Liouvillian space the right hand side of 
Eq.(S29) is a column vector. Therefore, our task is to 
solve a set of linear equations. At the end of the calcu- 
lation we fix the solution by projection out of the steady 
state by condition {{e\W)) = Ti-lTZJpst} = 0, accom- 
plished by premultiplication of projector Q|yV)). 



S. II. ADIABATIC ELIMINATION OF PHOTON 
DEGREES OF FREEDOM 



In the limit of strong photon relaxation, k ^ 7, g, the 
photon degrees of freedom decay so fast that the density 
matrix can be approximately factorized as 



p(t)^PDQD(t)(|0)(0|), 



(S30) 



where |0) is the vacuum state of photon in the transmis- 
sion line. Thus we can adiabatically eliminate the photon 
mode, and obtain the equation of motion for the reduced 
density matrix /5dqd(^, s) in the interaction picture' , 



PBQn{t, s) = riV{cl)pBQD{t, s) + TrT){cr)pBQn{t, s) + (7,. + 7ph)2?(cr )pnQB{t, s) + {s ~ l)J{(J )pBQn{t, s), (S31) 



S8 



where 7ph = ^g^/n. Note that the photon induced 
decay rate 7pij is associated with the spontaneous 
emission in this large n hmit, therefore photon ab- 
sorptions can be reflected by the jump superoperator 



J {a )/5DQD(i, s) = 7phO- PDap(i, s)ct+. In the ba- 
sis Pdqd = {po,Pg,Pge,Peg,Pe) : the matrix M{s) in 
Eq. (S15) is given by: 



M{s) 



( "4rj 

2Vi cos e + 2Vi 
-2Vi sine* 
-2Vi sin 6* 

V cose' + 2r; 



2Vr cos e - 2Vr 

— sin Q 
-r,-sin0 




— Fr sin d 
—Vr sin d 
-2r, - 27t 




—Vr sin Q 
—Vr sin Q 


-2r, - 27t 



2Vr COS 6 + 
Ajr + 4s7ph 



smf 
sin^ 



(S32) 



~r,.sin0 -r,.sin6' -2r,. cos - 2r^ - 47t / 



where ^ Jr + 7ph- 

To calculate the generating function G{s,t) and its 
derivatives, we take the Laplace transform of the gen- 
eralized density matrix, Eq. (S17), 



p(z,s) 



DQD 



(z-X(s))-lp(0,s)DQD. 



(S33) 



Since the long time behavior of the solution is determined 
by the residue of the generating function at the pole near 
z = 0, i.e., G{t,s) g(s)e^"* with g{l) = 1, we can 
expand the pole around s = 1: 



zo 



(S34) 



i>0 



and obtain, from Eq. (S22), the first two moments. 



■ Clt, 



s=l 



9^ 




[dg 




dg_ 






\ds 


s=l/ 


ds 



s=l 



(S35) 

(ci + 2c2) i, 
(S36) 



which give the mean and variance of the probability dis- 
tribution, respectively. In the asymptotic limit, t — > oo, 
all the information about the moments is included in the 
expansion coefficients c^. For instance, the Fano factor is 
given by"'^ 



(n) 



1 



2C2 
Cl 



(S37) 



To find the coefficients ci and C2, we consider the equa- 
tion 



det {zq\-M{s)) = 0, 



(S38) 



with Zo = ci(s - 1) + C2(s - 1)2 + O ((s - l)^). Then 
we can expand Eq. (S38) in powers of s and let the co- 
efficients for each power of s be zero. This procedure 
generates a set of equations with to arbitrary large i. 
We provide two examples below. 
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FIG. 5. (Color online) Dependence of photon noise Fano fac- 
tor on photon decay rate k at fixed bare coupling constant 
go = sr. Solid (dashed) lines represent results for numerical 
(analytical) calculations with method described in Sec. S.II. 
Other system parameters are Pj^^ = F, uo = 800F, T ~ 200F. 
Both plots indicate that the analytical results agree with nu- 
merical in the limit ^ go- 



First, considering Fj = F^ = F and — >■ tt with fixed 
coupling constant g, the case in which the two levels in 
the DQD are weakly overlapping, we obtain 



Cl 



C2 = -- 



7phr 

27tot +F' 
7^hF(7tot + 2F) 



F)' 



(27tot 

The Fano factor is given by 

p _ ^ 27ph(7tot + 2r) 
(27tot+F)' ' 



(S39) 
(S40) 

(S41) 



corresponding to the sub-Poissonian noise. When 6 = 
7r/2, a case when the DQD is tuned to its charge degen- 
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FIG. 6. (Color online) Dependence of photon noise Fano fac- 
tor on bare coupling constant go at fixed photon decay rate 
K = lOr. Solid (dashed) lines represent results for numerical 
(analytical) calculations with method described in Sec. S.II. 
Other system parameters are Fi.r = F, tJo = 800F, T = 200F. 
Both plots indicate that the analytical results agree with nu- 



merical in the limit 



, go- 



eracy, the solutions then read 

7phr(7tot + 2r) 



Cl 



C2 



67?ot + ll7totr 



4r2' 

2r) (47 



tot 



147t^ 



3i7totr2 



(S42) 
+ 20T^) 



(67tlt + ii7totr + 4r2)^ 



(S43) 



The Fano factor in this case is 



F = 1 



27ph (47L + 147?otr + 3l7totr^ 

(67?ot + ii7totr + 4r2)2 



2or^ 



(S44) 

again, giving the sub-Poissonian noise. In both cases, 
Fano factors are below 1 if 7ph is nonzero, indicating 
that it is the interaction between photons and electrons 
that gives rise to the sub-Poissonian statistics. 

As mentioned above, this analytical method is valid in 
the limit when k is large. We hereby make a comparison 
between analytical and numerical results, see FIG. 5 and 
6. In both plots, we do not consider the effects of dephas- 
ing on full counting statistics. The calculation indicate 
that the analytical method presented above agrees with 
numerical results in the limit k ^ go- 



S. III. CHARGE TRANSFER STATISTICS 

In this section we provide a quick review of the rela- 
tions for the electric current and current noise through 
the DQD using the master equation formalism. The cur- 
rent operator is defined as / = eF,.|i?)(i?| = eVrclcr- 
The dc current is given by the expectation value of / 
with respect to the steady state solution pst for the den- 
sity matrix, see Eq. (4) in the main part of the text: 



/ = eF,Tr{|i?) {R\pst}- 



(S45) 



The spectral density of the current fluctuations is de- 
fined by the relation 

/oc 
{{i{t)I{t + T)))e"^*dt, (S46a) 
-00 

{{iit)iit + r))) = (iit)iit + r)\ - I\ (S46b) 



The first term in Eq.(S46b) accounts for the concurrence 
of two electrons at times t and t + t 



{i{t)i{t + T))=l'gi',\t,T)+eIS{T), 



(S47) 



where the second order correlation function g^^ (r) is 
given by 

^(2)(^)^Tr{ctc.e^- [crP.cl]} ^ (S48) 



Tr 



1 4 Cr Pst I 



in terms of the steady state density matrix pst and the 
Liouvillian operator £, see Eq. (4) in the main text. The 
last term in Eq. (S47) represents counting the same elec- 
tron at t and t + t. 

Using Eqs. (S46) and (S48), we can write the current 
noise spectral function in the form: 



Siuj) = I' r {gf{T) - 1) e^-^dT + el. 



(S49) 



This result shows that the second order correlation func- 
tion .g^i (''') is related to the current noise in both fre- 
quency and time domains. In particular, the Fano factor 
Fci of the charge current that characterizes the low fre- 
quency limit of S{ui) is '^ 



el 



+ -/ (,g^f^(r)-l)dr. (S50) 



